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We study the dynamics of homogeneous, isotropic universes which are governed by the Eddington- 
inspired alternative theory of gravity which has a single extra parameter, Previous results showing 
singularity-avoiding behaviour for k > are found to be upheld in the case of domination by a 
perfect fluid with equation of state parameter lo > 0. The range — | < w < is found to lead to 
universes which experience unbounded expansion rate whilst still at a finite density. In the case 
K < the addition of spatial curvature is shown to lead to the possibility of oscillation between 
two finite densities. Domination by a scalar field with an exponential potential is found to also 
lead to singularity-avoiding behaviour when k > 0. Certain values of the parameters governing the 
potential lead to behaviour in which the expansion rate of the universe changes sign several times 
before transitioning to regular GR-like behaviour. 



I. INTRODUCTION 

In 1924 Eddington proposed a theory of gravity based 
solely on a connection field, without introducing a metric. 
A particular extension of Eddignton's theory including 
matter fields was recently considered in [T] (see also [5]). 
In a vacuum, this theory reproduces 'Eddington-inspired' 
(or just 'Eddington') gravity and is completely equiva- 
lent to General Relativity (GR), a fact which should be 
contrasted with other modified gravity theories in which 
vacuum deviations are expected (e.g. metric f{R) the- 
ories [3]). On the other hand deviations from GR are 
seen in regions of high density. Where one would ex- 
pect singularities such as the interior of black holes or 
the birth of the universe, in Eddington inspired gravity 
often these are avoided [U lU E] (although other singu- 
larities may appear in an astrophysical context The 
appearance of such behaviour in a purely classical theory 
such as Eddington gravity is interesting and exciting. 

The cosmological implications of this theory were ini- 
tially considered in [T], and studied in more depth in 
[5] and [7]; it was found that universes containing ordi- 
nary matter avoid singularites, with the exact behaviour 
depending on the sign of one of the parameters of the 
theory, and being either a regular bounce (from contrac- 
tion to expansion) or 'loitering' around a minimum size 
before transitioning to GR-like expansion. It has also re- 
cently been shown that tensor perturbations may render 
the regular behaviour of such cosmological models un- 
stable near the minimum scale (either in the asymptotic 
past or at the bounce [5]). 



In this paper we confirm those findings (in the addition 
of spatial curvature to the Friedmann-Robertson- Walker 
(FRW) metric), and present other interesting behaviour 
viz. universes which experience unbounded expansion 
rate at finite density akin to undergoing a second big 
bang/ 'cosmic hiccup', and closed universes which oscil- 
late between two finite values of the density. Finally we 
investigate the dynamics of universes in which the dom- 
inant contribution to the stress-energy tensor is due to 
a scalar field with an exponential potential, as has been 
done in the case of GR [5]. 

covers 



This paper is structured as follows: section I A 



some basic definitions and the governing field equations; 
section |I B| exhibits the Friedman equation and qualita- 
tively compares it to GR; sectionjlljexamines solutions to 
this with domination by one kind of perfect fluid; section 
|II B| demonstrates the 'loitering' behaviour which has pre- 
viously been observed; section [II C| investigates the 'cos- 
mic hiccup' mentioned above; section II D exhibits the 
oscillating universes; section [Tll] cove rs the case of domi- 
nation by a scalar field; section III A| investigates the 'loi- 
tering' behaviour which has been observed for the case 
of regular matter fields and which shows up again in the 



case of scalar fields; section III B exhibits some universes 
which show qualitatively new behaviour - expansion, fol- 
lowed by contraction, followed by expansion. In section 
HVl we conclude. 



A. Definitions and Field Equations 

The action governing Eddington gravity is [I] 
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(throughout this paper we use reduced Planck units, c = 
SttG = 1) where R{T) denotes the symmetric part of 
the Ricci tensor, and A is a dimensionless constant which 
must satisfy A = 1 + kA if GR is to be recovered in the 
hmit of kR ^ g (though the cosmological constant will 
be assumed negligible in the following). Matter can, of 
course, be coupled to the connection as well as the metric, 
however for simplicity here it is assumed to only couple 
to g. Two field equations are derived from this action: 

qt^y = 9t,y + nRt,i>{q), (2) 

v^?''" - ^Mg^"" = (3) 

in which q is an auxiliary metric which is compatible with 
the connection i.e. T^^ = ^q'^^idpg^^ + d^qp^ - d^qp^), 
and q^''' = [q-^f . 

The action (Tj), and its equations of motion (2|3), can 
also be derived from the following bigravity theory, 



I[9. q] 



with 



V^(q'"'i?^.(g)-2A) (4) 
- Ki^q^'^g^,, - 2^g) + ^L,n{g) 



A=- (5) 



Varying Q with respect to g^i, and g^jy one obtains a set 
of equations fully equivalent to (2j3|. This form of the ac- 
tion is particularly interesting because it shows a remark- 
able close relation with the recently discovered family of 
unitary massive gravity theories ^lOj . Indeed, massive 
gravities are built as bigravity theories where the poten- 
tial is a linear combination of the symmetric polynomials 
of the eigenvalues of the matrix 



= Vq ^9- 



(6) 



For this class of potentials the theory has no Boulware- 
Deser instability and, in particular, the massive sector 
carries 5 degrees of freedom instead of 6. Besides the 
cosmological constants (which also appear in Q), one 
such symmetric polynomials is 

(Tr(7^) - (Tr7)^) = V^iq^'g,, - {Trj)') (7) 

The first term is precisely the interaction between q^^, 
and 3^1/ appearing in Q, resulting in Eddington theory 
coupled to matter. The second piece (more difficult to 
compute because the square root has to be taken before 
the trace) is not part of the Eddington theory. We find it 
intriguing that Eddington's theory coupled to matter is 
so close to the recently discovered stable massive gravity 
theories. We shall exploit this connection elsewhere. 



B. Friedman Equation 

The FRW metric is 

ds^ = -dt^ + a{tf ( + r\d9^ + sin20d02) 



and we take the auxiliary metric to be 

q^^dx^'dx" = - U{t)df + a{tfV{t) 
dr^ 



1 - fc,r2 



r^(dr -l-sin-'e'c 



(9) 



The energy momentum tensor is 




Too ^'^Pi + P4' and T^^ 

(10) 

where {i} indicate the regular matter, and (fi the scalar 
field. Finally 

P^^l^' + V{^) and P^ = ^^'-V{cl,). (11) 

The Eddington gravity field equations enforce kg = k 
and 



C/=W^V^^ and V=^{1 + >,pt){1-kPt), 
y l-f npT 

(12) 

where Pt = J2i Pi + P'P + ^ and Pt = I]* Pi + P<j, - -A- 

The energy-momentum conservation equations are as 
in GR, so 



(l3 + 3H(l3 + V'{(j)) =0 
p^ + 3Hp,{l + w,) =0. 

In the general case 



H 




B 



with 

G- 

F 



= 1 



--4 



duk 



(13) 
(14) 



(15) 



(16) 



Pt 



B = 



1 

2'^1 - kPt' 



Wi ~ K{wiPT + Pt)) Pi(l + Wi) 
' 4{1 + kpt){1 - kPt) (17) 

(18) 



(8) 



One key thing to note is that the scalar field does not 
simply enter in the same way as other forms of matter (as 
is the case in GR), which leads to the interesting result 
that although there are two solutions for H, they are not 
simply opposite signs, as in GR and Eddington gravity 
without a scalar field. 

If one assumes that sign(V^) — — sign(i?), i.e. as the 
universe increases in size, the potential energy of the 
scalar field decreases, then in fact the values for H are 
simply of different signs. However whilst this assumption 
may hold for the low-density GR-like stage of evolution. 
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it fails in the high density case with some interesting con- have to have in order to evolve as curvature) . 



sequences as laid out in section IIIB 



Another departure from GR is that spatial curvature 
cannot be considered a fluid, since it does not enter in 
the same way as a fluid with w = (which it would 



Finally, in the low density and curvature limit 
(pi, A, fc/a^ <C K^^) GR is recovered to leading 
order: 



H =± 



± 



Pt 
3 



-E 

4 



12 I 3 



15p2 9 17 2 

— Pj. - -ptPt - 



-^(Pt + Pt, 



(19) 



II. DYNAMICS OF A UNIVERSE WITHOUT A 
SCALAR FIELD 

We first consider the dynamics of universes which do 
not posses scalar fields, have a negligible cosmological 
constant and, for simplicity, only one type of perfect fluid. 
The Friedman equation in this case is 



G 
8 

3k 



(20) 



(1 + ■iw)p - 2 + 2^J(l + p){l-wpY 

- J(l-^p-) 

(l + p){l-wpf 

(4 + (1 - w){l - 3w)p + 2w{l + 3w))^2)2 



(21) 



1 (so that the universe con- 



where p = Kp and k = Kk. 

Note that setting w — - 
tains only a cosmological constant) reduces eqn. (21) to 
exactly that which would be expected in GR, as it should 
be since Eddington gravity exactly reproduces GR in a 
vacuum (with a cosmological constant). 



A. Behaviour of the Friedman equation 

The square root in the numerator restricts the domain 
of and the result is shown in table |ll 



TABLE I: Domain restrictions placed on H by the square 
root 



w 


Valid 


P 


range 


w > 


-1 < 


P 


< i 

— W 


w = 


-1 < 


P 




-1 < w < 


P<- 


or 


-l<-p 


w <-l 


P<'l 


or 


-<P 



TABLE 11: Sign of p in various regions of parameter space, 
with no spatial curvature, pc* satisfies G{pG,) — H^{pG*) = 
0. {w — —1 is special but also has p > 0.) 



P 



sign(K) sign(p) 



any < p 

-1 <w -1 < p < 
w < pa* < p < ^ 
P < PG* 
W < ~1 pG* < p < 
i<P<pG, 
p< -1 



+ 



+ 

+ 
+ 

+ 



The requirement > fixes the sign of p and table 
[TT| shows the results (for fc = 0) , from which we see that 
that the energy density is positive in all but two regions 
of parameter space. This is an interesting departure from 
GR which, in the spatially flat case, always requires non- 
negative energy densities. 



As pointed out in 7], a final interesting quality of (21 ) 
to note is that for certain ranges of w, — > oo for finite 
p. Section [lI C| investigates a universe which exhibits such 
behaviour. 

Looking at the discriminant of the quadratic in the de- 
nominator of (21), we see that the discontinuities do not 
exist for 0.02 ^(7 - 4^3) <w< 1(7 + 4^3) w 4.64. 
Outside of this region we can determine graphically, given 
the constraints of table |Tj in which regions the disconti- 
nuities will appear (see appendix [A|) , and the results are 
in table Unl 

Points at which ~ and ^ will induce an 
expanding universe to collapse (and vice versa) since the 
universe cannot evolve further on this trajectory (because 
> 0) and the stationary point is unstable since only 
if^, and not the sign of H, is fixed by the physics. Such 
'bounce' points exist in GR (if fc > 0), and in Eddington 
gravity as well, i.e. the points p = — 1 and p = pc*- 
We are not including p = here because, although 
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TABLE III: Regions in which discontinuities of appear and are not outside of the allowed range of p. 



w<— 1 w = —1 —1 < w < — \ w = —\ 


-I <w <0 0<w < 4.64.. 


4.64... < w 


pa* < p <0 - p < pG* 


< p 


-1 < p < 



behaves as required, starting from a non-zero density this 
point ca nno t be reached in a finite time. 

From (j2lj) we see that at p = w^^, iJ^ = Q and ^ = 
0. This novel type of stationary point leads to the 'loiter- 
ing' behaviour remarked upon in pQ |S] , and [7] , in which 
the universe asymptotically approaches (or recedes from) 
a maximum density/minimum size. From table [ll] we see 
that for a positive energy density, this type of universe 
will only exist for w > and k > 0. 



B. 



The case k > 0, w > 



As seen in figure [T] this type of universe exhibits the 
'loitering' behaviour mentioned in the previous section. 
We also see that to a good approximation the evolu- 
tion can be split into two phases: the loitering phase in 
which deviations from GR are important, and the GR- 
like phase in which the dynamics of the expansion are 
approximately as given by GR. 



Calculating the particle horizon for such a universe we 
find 



rp(t) = lim 



dt' 



t-^-oo J^.^ a{t') t 



lim 



(26) 



which diverges. Again this is similar to GR under domi- 
nation by a cosmological constant. 

In [3] the authors note that Eddington gravity can be 
characterised as making the gravitational force repulsive 
at short distances (or high densities), which would agree 
with the behaviour which we see here. That is, the max- 
imum density is reached as the repulsive effects of Ed- 
dington gravity oppose the traditional attractive effects 
of the matter density; which repulsive effects take a form 
similar to a cosmological constant with A = 8k^^. 



2. Behaviour of the auxiliary metric 



1. Loitering phase 

Close to the maximum density p = ps = w 
involving k is much smaller than the other terms in the 



numerator of (21 1, and so we can neglect (small) spatial 
curvature during the loitering phase. Expanding about 
the maximum density, or equivalently the minimum scale 
factor as, to highest order the Friedman equation and 
deceleration parameter become 



Qde 



wSp 



3k V3(1 + w) 
w6p ^ ^ 



3(1 



8 fSa 



(22) 
(23) 



(where 5p 



5a 



Pb and 5a = a — clb)- Using H = ^ = 



( 22 ) can be solved to get 



a(t) 

Qb 



qde 



= 1 



CXp I \j 7^it-to) 



-exp 




to) 



(24) 



(25) 



Thus we see that the universe undergoes (accelerating) 
exponential expansion away from the minimum scale fac- 
tor during this phase. The form of a(t) is reminiscent of 
that observed in GR under domination by a cosmological 
constant, specifically with A = 8k~^. 



From eqn. (j9| we see that the auxiliary metric (7^,^ is 
identical to the metric g^^ except that the temporal and 
spatial components of q are modulated by the functions 
U and V respectively. Alternatively, by factoring out U, 
it can be considered to be related to an FRW metric with 
modified scale factor ag = a^^' overall conformal 

transformation U. Thus U and V, or U and Uq explain 
how the space described by the auxiliary metric behaves. 

Figure [T] shows the behaviour of U, and V. As ex- 
pected, in the GR regime U, V ^ 1. On the other hand, 
close to the maximum density U, V ^ 0, and the space 
described by q becomes vanishingly small. Interestingly, 
V reaches a maximum value just before the GR regime 
begins. 

Expanding close to the maximum density we find 




3(1 + w) 



(27) 
(28) 
(29) 



And so we see that, moving towards the maximum 
density, the temporal size decreases more rapidly than 
the spatial, with the result that Og actually goes infinite. 
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FIG. 1: Left: Expansion rate of the universe against density (normalised by the maximum density pa — (kw) for k — I, 
w = 1/3, fc = 0. 

Right Top: Scale factor (normalised by the minimum scale factor as) against time (in reduced Planck units), for the full 
numerical solution to the Friedman equation, and approximations in two phases: loitering, as in section |1I B 1| and GR-like 
low-density limit. 

Right Bottom: Behaviour of the parameters describing the auxiliary metric (for the same universe). 



3. The case w = 

This case exhibits the same loitering behaviour at early 
times, but must be treated separately, since the density 
is unbounded; expanding about a = 0, to leading order 
we find: 



3k 

-1, 



(30) 
(31) 



which is identical to GR with a cosmological constant 
A = 8k~^. Therefore the particle horizon will diverge as 
in the case w > 0. 

The behaviour of U and is the same as in the w > 
case, however in contrast, V diverges as t — > —oo. See 



appendix B 1 for the explicit functional forms. 



C. The case k>0, -|<it;<0 

Firstly it should be noted that new physics (beyond 
Eddington gravity) is required in order to have a fluid 
with w < dominating in the early universe; this prob- 
lem is not considered here. 

As shown in table [ITU such a universe the expansion 
rate diverges at a finite density given by the positive root 
of the denominator of eqn. (21 1 



(1 - ii;)(l - 3w) + (1 + w)\/l -42u; + 9w2 
4|w|(l + 3w) ■ 



(32) 



Expanding around this point, to leading order we find 

-2 



OK \ a* 



(33) 



with 



(1 + 3w)p* - 2 + 2v/(l + p*)(l- wp*)^ 
6~k _ 1 1 



(3p,(l + w;))2 
(1 + p*)(l - wp*)2 



((1 - w)(l - 3w) + 4w(l + 3u;)p,)2 ■ 
This can be solved as in section [IIB l| to get 



a(i) 
a* 



1± 



l2U\j—\t-U\. 



(34) 



(35) 



Thus at this critical point the universe in some sense 
undergoes a second big bang, with expansion similar to 
that of GR with domination by radiation (a cx ^/i). Fig- 
ure [2] demonstrates this universe. 



Expanding eqn. (21) for high densities we find 



\w\ 2 



3 (1 + 3^)2"^- 

The expansion in the early history of this universe is also 

3 

GR-like, albeit with a modified density p = jj^f^^^P, 
and the spatial curvature can be neglected. The scale 
factor as a function of time is then 



(37) 
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FIG. 2; Left: Expansion rate of the universe against density (normalised by the critical density p,), for n = 1, w — —1/6, 
fc = 0. 

Right Top: Scale factor (normalised by the critical scale factor a*) against time (in reduced Planck units), for the full numerical 

solution to the Friedman equation, and an approximation close to the critical density. 

Right Bottom: Behaviour of the parameters describing the auxiliary metric (for the same universe). 



In order to better understand the dynamics of the uni- 
verse during this second big bang/ 'cosmic hiccup' we 
should consider the horizon between two times t_ < < 
t, close to i*: 



high densities more quickly than — > oo, as indeed it 
should if the high density behaviour observed in section 
III B 31 is to be recovered when w — 0. 



r(t_,t+) 



t+ - t_ 



(38) 



Its well-behaved nature should not be surprising, how- 
ever, since the nature of the expansion during the 'hic- 
cup' has already been noted. 



1. Length of first life 

An interesting quantity to consider is t*, the time be- 
tween the birth of the universe and the 'hiccup'. Naively, 
since — ^ oo for w — > 0, — |, it might be expected that 

— >■ in these limits. 



We have 



da 



aH{a) 



dp 



pH{p) ■ 



(39) 



Whilst this cannot be evaluated analytically in the gen- 
eral case, certain approximations can be made. For 
w — )■ 



p* cx 



l+3u) 



, and H is given by eqn. (36) and 



so the behaviour of eqn. (39) can be approximated by 

« C = "^gr^^ (1 + ^^f" ^ 0: as ex- 

pected. 

However for w — > 0, p* 



haves: t* 



CP, 
J oo 



dp 



1 

2\w\ 



and so eqn. (|39|) be- 

oo. 



1-1/4 



-loo \wp/'^p^/^ ~ |^|3 

which defies the naive expectations. This occurs because 
the expansion rate of the universe heads towards zero at 



2. Behaviour of the auxiliary metric 

In the high density limit the overall size of the space 
described by q diverges, but the combined effect of U and 
V is that Uq remains proportional to the regular scale 
factor. 

At the 'hiccup' U and V as functions of p or a are 
smooth and well behaved, but as functions of t they of 
course exhibit a 'kink' just as a{t) does (see figure [2]). 

Explicit functional forms in these two limits can be 



found in appendix B 2 



D. Oscillating universes 

In [TJ [S] and [7] it is noted that for k < 0, w > 
the universe undergoes a 'regular' bounce (as described 
below) in which, in a finite time, a collapsing universe 
reaches a maximum density and then expands. In GR, 
universes with positive spatial curvature can undergo a 
bounce in the opposite direction. 

The addition of spatial curvature to Eddington cosmol- 
ogy thus leads to the interesting prospect of oscillatory 
behaviour. Figure [3] demonstrates such an oscillatory 
universe. 



7 



0.08 




0.4 0.6 
P/PB 



ml.5 

C3 



1 

1.5 

1 

0.5 




1 1 1 1 1 1 1 1 


\ 1 
I 


1 \\ 1 




- / Exact Soln. 

-A High Density Approx. 

GR Approx . 


\v 




- 


J \|^^^^'^' ^ ' ' 


— 1 


1 y^l^* 


H — — : 


















i — u 
/ V 


■ 1 V 


1 1 1 


. 1 - 







10 



15 



t 



FIG. 3: Left: Expansion rate of the universe against density (normalised by the maximum density pB ~ — k for k — — 1, 
w = 1/3, k/a% = 1/10. 

Right Top: Scale factor (normalised by the minimum scale factor as) against time (in reduced Planck units), for the full 
numerical solution to the Friedman equation, and approximations in two phases: close to the maximum density, and GR-like 
low-density limit. 

Right Bottom: Behaviour of the parameters describing the auxiliary metric (for the same universe). 



1. Amplitude of the oscillations 

The maximum density, ps = is independent 

of w and fc, and is caused by the (1 + p) factor in the 
numerator of eqn. (21 1. The minimum density, on the 



other hand, does depend on w and k (and k), and is the 
solution to G = 0, which explicitly is 



(1 + iw)p - 2 + 2^ {I + p){l - wpf 

- 6fca"2(l - wp) = 0. (40) 

Unfortunately this cannot be solved analytically in the 
general case, however if we assume that the universe has 
expanded sufficiently before expansion halts then we can 
solve this in the GR regime to get 



3fc 

7,2" 



Pmin 
PB 



3k 



3(1 + ^) 
1 + 3™ 



(41) 



From figure H we see that this is indeed a good approx- 
imation for fc/a^ ^0.2. 

Close to the maximum density we have 



31^1 

which is solved to give 



1 + 



2k\ Sa 



ait) 

Qb 



1 



2 

3M 



1 



2k 

7X 



\t~ti 



(42) 



(43) 



Note that this is independent of the t ype of m atter which 
fills the universe (just as in section IIBll, and so we 



understand the expansion is being driven by k. 



Bounds can be placed on k by examining the Friedman 
equation in these two limits (since the behaviour of eqn. 
(211 as a function of p and k is smooth) and requiring 
0. The GR-like regime gives a lower bound on k 
(42) gives an upper bound, which combine 



whilst eqn 
to give 



< 



k 



< 



(44) 



2. Period of the oscillations 

In principle the period of the oscillations can be calcu- 
lated from 



Tosc = 2 dt^2 
Jo 



da 



aH{a) 



(45) 



In practice this integral cannot be performed analytically 
(or good enough analytical approximations made) and so 
must be evaluated numerically, as shown in figure [4] 

We see that Tosc{w = |) < Tosc{w = 0), which makes 
sense since the amplitude of the oscillations for the for- 
mer case is always smaller than in the latter. A lso note 
that Tosc for k/a% — >■ 1/2. From eqn. (42 1 we see 



that the expansion rate goes to zero at this point, and 
so from section II CI we might expect that Tosc 



oo. 



however in this case the extra size by which the universe 
must expand also heads to zero (and clearly does so more 
quickly than the expansion rate). 
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FIG. 4: Top: Amplitude (minimum density, normalised by 
maximum density) of the oscillations described in section pl D| 
against spatial curvature of the universe, for k = — 1 and vari- 
ous values for the equation of state parameter. Full numerical 
solutions, as well as approximations in the GR regime from 
eqn. ((411. 



3. Behaviour of the auxiliary metric 



Close to the maximum density the spatial size de- 
creases to zero as the universe reaches is minimum size, 
whilst the temporal component of the auxiliary metric 
diverges. See appendix |B 3| for approximate functional 
forms for these. 



Bottom: Period of the oscillations (in reduced Planck units) 
for the same universes. 



III. DYNAMICS OF A UNIVERSE WITH A 
SCALAR FIELD 



We now consider the cosmology when a scalar field is 
introduced; the ordinary matter content, the cosmolog- 
ical constant, and the spatial curvature are considered 
negUgible. Written in terms of 4> and V{(j)) the FRW 
equation is 



H = 



± [1 - -<j)^ + kV 



\ 













1 









1 + W 1 + + 1 - -02 + 



1 -f 2kV - 



(46) 



The positive sign is taken, as this gives a universe 
which is expanding in the GR regime. The conservation 
equation governing the scalar field is 



(47) 



We restrict ourselves to an exponential potential for the 
scalar field 



V^(</)) = yoexp(-A</)) 



(48) 



with Vo, A > 0. Unless explicitly included, the parameter 
describing deviations from GR (k) is set to one. 



A. Loitering behaviour 



For K > numerical solutions of eqn. ( 46 1 show qual- 



itatively similar behaviour to the case without scalar 
modes, viz. 'loitering' of the scale factor around a mini- 
mum value followed by a transition to GR-like power-law 
expansion. See figure [Sj top left panel. 



The solution for 0(t) also shows two regimes with a 
transition between them (figure [sj bottom right panel) . 
The GR regime shows the expected logarithmic time 
dependence [9,, whilst the Eddington/loitering regime 
shows cx —t. 



1. Initial conditions and free parameters 

It is appropriate first to briefly discuss the initial con- 
ditions used to produce the solutions plotted above as 
well as the effect on these that is had by altering the free 
parameters of the system. 

As can be seen from figure [5] the chosen conditions are 



a{t = 0) = 1; 0(0) = 0; 0(0) = 0. 



(49) 



The first is uncontroversial since the scale factor can al- 
ways be rescaled by a constant. The absolute value of 
the scalar field only appears in the form for its poten- 
tial energy, and thus adjusting the field by a constant. 
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FIG. 5: Scale factor (top left), scalar field (bottom left), energy density (top right) and equation of state parameter (w = P/p) 
(bottom right) as a function of time (in reduced Planck units) for various values of the free paxameters. 



which is the action of the second i.e., only has the ef- 
fect of changing Vq. The theory is obviously invariant 
under a constant time shift and the third i.e. just fixes 
this. The reader may be worried that using = as 
a i.e. eUminates from consideration any solutions which 
show monotonically increasing field values, however such 
solutions would not exhibit the maximum density and 
loitering behaviour (since if (p is unbounded from below 
the density will also be unbounded, and if (j) tends to 
some constant value, the pressure would be negative and 
hence not satisfy P = k~^). 

All that's left is to describe the potential, which leaves 
two free parameters: Vq and A. Prom the figures thus far 
presented we can roughly evaluate the effect of making 
the potential steeper (increasing A) or deeper (increasing 
Vo). 

The late time (GR-like) behaviour is primarily gov- 
erned by the value of A, whereas the early time behaviour 
is infiuenced by Vq. The reason for the latter is that the 
energy density of the field is already fixed as t — >■ — oo, 
and picking a value for Vq fixes it at t = (or wher- 
ever the transitionary period is located), and thus the 
behaviour in between these two points will be greatly in- 
fluenced by the size of the change in the energy density. 

Furthermore one can roughly say that steeper and 
deeper potentials lead to quicker transitions from the loi- 



tering to the GR phase. 



2. Loitering regime 

Let us investigate this loitering behaviour in more de- 
tail. As in the case with ordinary matter the loitering 
regime corresponds to P — >■ k"^ (as t — >■ — oo); we also 
see that — >■ oo and hence V{(j)) — >■ 0. Thus one has 

P and (f> ~'\f^i + const. 
Evaluating the Friedman equation in this regime, one 
finds 

HiP = .-^)='^^±^='V2nV'. (50) 

Therefore the equation of motion for the scalar field re- 
duces to ^ + 3V' = 0. The solution (setting k = 1) is 

^(t) = -y/2{t - to) + ^ 21n(12yo + e^^(*-*°^) - 4 In 2 

(51) 

a{t) w amin exp ^— ^ exp \/2A(i - to)^ . (52) 
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The dynamic nature of the scalar field adjusts the loi- 
tering behaviour of the scale factor from the simple ex- 
ponential behaviour seen in the case of regular matter to 
the more extreme double exponential. 



3. Transitionary Regime 

Whereas in the case of ordinary matter fields the evolu- 
tion is well described without having to denote a separate 
transitionary regime (between the high-density loitering 
behaviour and the low-density GR-like regime), that is 
not the case here. 

The transition from linearly decreasing field values to 
logarithmically increase is naturally centred on (j) ~ 0. 
Hence during the transition P « —p: the scalar field has 
the equation of state of a cosmological constant. This 
is readily observed in plots of pit) which show that the 
density is stationary around the transition. See figure [5] 
top right panel. 

The bottom right panel of this figure also confirm this, 
as well as the above comment that P « p in the loitering 
regime [w — )• 1 as t — > — oo, regardless of the values of 
the free parameters). 



4- Density Scaling 

Finally it is interesting to look at the plot of In p 
against In a (figure [6| . This shows that during the loi- 
tering phase the density scales as p oc a~^, which agrees 
with the fact that w 1 - the scalar field is, to use 
the terminology used in [9 , undergoing kination (kinetic 
energy domination). This is followed by the expected 
cosmological constant-like transitionary phase in which 
the density is stationary. The field energy density shows 
the expected scaling in the GR phase (p oc if A < 
•\/6 and p oc otherwise). 

Figure 1 from |9] shows the evolution of the energy 
density of a scalar field in a GR-universe populated by 
the aforementioned scalar field (with A = 4) as well as 
radiation and cold matter. The scalar field is initially 
dominant and scales accordingly (p cx a~^), becomes sub- 
dominant and scales like a constant, before finally track- 
ing (i.e. scaling as) the cold matter (p cx a~^) which is 
then dominant. 

This bears a remarkable resemblance to figure[6j except 
for two major differences: 

1. the behaviour in the Eddington case is reversed. 
The 'natural' scaling (i.e. that which one would 
expect in GR based on the value of A) is preceded 
by a period in which the scaling is fixed (by 
something other than A); 



2. in the Eddington case such behaviour is achieved 
without another fiuid for the field to track. 




Vo=l/4',A^2' 
Vo=l/4, 7^3 ■ 
y 0=1/2, 7^2 
Vo=l/2, >^3 
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FIG. 6: Energy density as a function of scale factor for various 
values of the free parameters. 



B. Expand-contract-expand behaviour 



As noted in section |IB[ the two solutions for the Hub- 
ble parameter are not in general simply of the same mag- 
nitude but with opposite sign. This leads to the possi- 
bility that there are universes which at some point are 
stationary, and then instead of contracting and retracing 
their evolutionary history in reverse, they contract but 
trace a new history. 11 

In fact since the solutions being investigated are ex- 
panding both in their early and late histories, they would 
undergo expansion, followed by contraction, followed by 
expansion. 

From figure [7] we see that such behaviour can be 
achieved both by making the potential steeper and by 
making it deeper. The latter can easily be understood 
since Vq, which controls the depth of the potential, also 
controls the energy density of the scalar field at its tran- 
sitionary point. Thus choosing a sufficiently large Vq will 
ensure that the density at this point is larger than the 
density in the distant past (loitering phase), and hence 
the universe must undergo contraction (see figure [7j bot- 
tom panel). 

That steeper potentials induce this behaviour can be 
understood mathematically from the fact that the 'B' 
term in the Friedman equation (eqns. 15 18) is propor- 
tional to V' , and hence to A. Physically, perhaps one 
can imagine the cause to be that a potential which varies 
quickly enough with (j) leads the scalar field to 'overshoot' 
the value which would ensure energy density monotoni- 
cally decreases with time; this larger density would then 
be associated with a contraction. 

It is interesting to note that plots of the field value, or 
of the equation of state parameter do not obviously show 
any affect of the period of contraction or differ qualita- 
tively from the corresponding plots for free parameter 
values which show the ordinary behaviour (figure [s]). 
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t 

FIG. 7: Hubble expansion parameter (top), scale factor (mid- 
dle), and energy density (bottom) as a function of time (in 
reduced Planck units) for various values of the free parame- 
ters, showing universes which expand-contract-expand. 



domination by a perfect fluid is upheld for w > and that 
this result is not affected by the presence of spatial curva- 
ture. This behaviour shows remarkable similarities to the 
behaviour in GR with a cosmological constant A = 8k~^. 
Similar behaviour is also seen in the case of domination 
by a scalar field with an exponential potential. The scal- 
ing of the density with the scale factor in such universes is 
reminiscent of that observed in a GR universe with such 
a scalar field and (initially more slowly scaling) ordinary 
matter. The similarities to GR in each of these cases de- 
serve further investigation which will shed light on the 
high density behaviour of Eddington gravity. 

For both perfect fluids and scalar fields we have also 
demonstrated novel behaviour unlike that seen in GR. 
In the former case we have the 'cosmic hiccup' which is 
observed for k > 0, — | < w < 0, as well as the oscillatory 
universes for k < 0, w > 0, and positive spatial curvature. 
In the latter case we have shown that especially deep or 
steep scalar field potentials lead the size of the universe 
oscillating whilst transitioning between the loitering and 
GR-like phases. 

There appears no initial reason why the singularity 
avoiding behaviour observed for k < cannot also exist 
for a scalar field dominated universe. As explained above, 
symmetry around the bounce would be permitted in this 
case (0 = at the bounce). Furthermore the interplay 
between the scalar field and regular matter fields could 
lead to interesting behaviour of the scalar field, as is the 
case in GR [51 and remains to be studied in detail. 



Unfortunately we have yet to determine the critical 
values of Vq and A which guard the transition into such 
behaviour. 



IV. CONCLUSION 

In this paper we have shown that the singularity- 
avoiding 'loitering' behaviour observed for k > and 
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Appendix A: Regions in which discontinuities appear 

Figure [8] shows, for a universe dominated by a single 
perfect fluid, the locations in (p, w) parameter space of 
discontinuities in H^. There are no discontinuities at all 
for 0.02 « 5(7 - 4a/3) < u; < ^(7 + 4v^) « 4.64, and 
outside of this range they are given by 



p* 



(1 - w)(l - Sw) ± (1 + - 42u; + 9^2 

4w(l + 3w) 



(Al) 

Also shown are the regions which are forbidden, due 
to the square root in G (and the requirement that is 
real). Combining these data we arrive at table III 



Appendix B: Approximate functional forms for the 
behaviour of the auxiliary metric 

1. The case k > 0, w = 



Where a is given by (37) and thus 



The functional forms close to the critical point are 



V = v. 



1 - 



1 - 



3p*(l + w) f 1 



3w \ Sa 



1 + p* 1 — wp^ J a* 

3p^,{l + w) /I w \ Sa 

2 \ 1 + 1 — wp^ J a* 

3p*(l + w) / 1 w 



I + p* 1 - wp+ 



and — can be read off from ( 35 1 



(B5) 



(B6) 
Sa 
a* 
(B7) 



Defining p(a = ao) ~ po, during the loitering phase we 
find 



U 



1 fSaV 



Po \ao/ 



1 / 3 

: exp 



2 V 3k 



it -to)] (Bl) 



'Po 



= VPoGXp|--Y-(t-to) 



aq = aoVPo 



aoVPo exp 



2 V 3k 



(B2) 

(t - to)^ 
(B3) 



2. The case k>0, -|<«;<0 

In the high density limit, for the early history of this 
universe, we have 



3. The case k < 0, u; > 



Close to the maximum density we have 



U = (1 + w), 



V ={l + w) 




t-tf 



(B8) 



+ + ^ ]\t-tB\ 



3<5a=^(l 



2k 



\t~tB\ 



(B9) 
(BIO) 



U^\w\^p, V=\w\^p, aq = a\w\~ 



(B4) 




FIG. 8: Hatched regions are those forbidden due to the square root in G; sohd lines indicate the solutions to F = and so 
represent the locations of discontinuities in H^ ; dashed lines indicate the asymptotes to which the solid lines tend. 



